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1. Introduction 

The theoretical description of decay processes which include light mesons strongly relies on a 
non-perturbative and model independent determination of the relevant parameters. The experimen- 
tally most accurately known decay constant is the pseudo-scalar decay constant, /ps, which can be 
determined from the leptonic decay of charged pions or kaons. Attempts to compute this quantity 
from first principles on the lattice is an important validation of the lattice QCD methodology. It 
allows in particular to test control on the required extrapolations to the infinite volume, continuum 
and finally chiral limit. For phenomenology much more interesting is a lattice determination of the 
coupling of the light vector mesons to the vector current (/v) and the tensor current (/y )■ While 
the constant /v can in principle be extracted from experimental studies of, e.g., leptonic decays like 
-^V^V-i, the coupling /y can only be determined theoretically. It has been pointed out by the 
authors of [[T]] that the form factors describing the semi-leptonic B-decay B —>■ plv strongly depend 
on these couplings. A precise calculation of these couplings is hence an important step towards a 
more precise determination of one of the less well known CKM parameters, |Vub|- 

2. Lattice Details 

The configurations with N( = 2 flavours of non-perturbatively improved Wilson fermions and 
Wilson glue which have been generated by the QCDSF and UKQCD collaborations in recent 
years now cover a larger range of quark masses and several lattice spacings. The simulations 
are performed for several values of the gauge coupling jS with the lattice spacing varying between 
0.07 — 0.11 fm. For each lattice spacing configurations are generated for 3-4 sea quark masses in 
the range of 590 — 1 170 MeV. For additional leverage on the quark mass dependence we performed 
simulations at a large number of valence quark masses which vary between 470 and 1200 MeV. 
Also the volumes are varied to explore finite size effects. Finally, a non-perturbative determination 
of the renormalisation constants is performed to further reduce systematic uncertainties. 

Correlation functions for pseudo-scalar and vector mesons ai^e calculated on configurations 
taken at a distance of 5-10 trajectories using 8-4 different locations of the fermion source. We use 
binning to obtain an effective distance of 20 trajectories. The size of the bins has little effect on 
the eiTor, which indicates auto-coiTclations ai^e small. To improve overlap with the ground state we 
employ Jacobi smearing. 

In order to scale our results at different values of j8 and K"sea we use the Sommer-scale ro, 
which can be determined with good precision on the lattice. When translating this into more com- 
mon physical units, we face the problem that the static quark potential is difficult to determine 
experimentally. We therefore choose to use lattice results for rowiN (where win is the nucleon mass) 
to obtain ro = 0.467 fm. We typically will use r^/a at the respective values of j3 and Jfsea- Only 
when stated explicitly we will use r^/a extrapolated to the chiral limit [^. 

To fully eUminate 0{a) discretisation errors, it is not sufficient to improve the action, the 
operators also need to be improved. To determine the renormalised pseudo-scalar decay constant 
we therefore have to calculate 




(2.1) 
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Figure 1: Comparison of ca determined from perturbation theory (dotted line, [g]), boosted perturbation 
theory (dashed Une), and non-perturbative methods (sohd line, |W|)- 



where Za is the axial vector cuiTent renormalisation factor and bp^ a coefficient describing its depen- 
dency on the bare sea-quark mass anif^^^sa = 5(1 /'^sea — 1/ '^fea)- '^a is the improvement coefficient 
which needs to be chosen such that 0{a) discretisation errors cancel. 

The unimproved bare decay constant f^^ and the improvement term are computed from the 
matrix elements mpsfp^ = (0|A4|PS) and mpsafp^^ = {0\ad4P\PS). While the renormalisation 
constant Za [||, (we use [||]) and the the critical value of the hopping parameter k^^^ [^] have 
been determined non-perturbatively, the coefficient Z^a is only known perturbatively [^. As the 
expansion in the bare gauge coupling = 6/[5 is known to have bad convergence properties, 
we use the boosted couphng ^^■^(ju) at the scale jj. = I /a. This couphng is computed from the 
renormalisation group equations with n/A^^ = {vq/o) / {r^h!^^), taking vq/q in the chiral limit 
and vqA!^^ = 0.617 (40) (21) [||]. To assess the differences between perturbation theory, boosted 
perturbation theory and non-perturbative methods, we compare the results for ca- As can been seen 
from Fig. H the values from boosted perturbation theory aie much closer to the non-perturbative 
results. However, there is still a significant discrepancy between the two results. 

Similarly, we define the renormalised and improved vector decay constant /y: 
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(2.2) 



where e{X)im\ fy^ 



(0|y,-|V,A) 2inAe{l)im\af'^'^ 



/V ~ " '"•/ ^\'^ ji'i^y^jM — (0|a547^'4|V, A) with e(A),- being the po- 
larization vector. The renormalisation constant Zy and the coefficient b\ have been computed 
non-perturbatively [0], but we have to rely on a perturbative calculation of cy [^. 

Finally, we define the renormalised and improved constant for the coupling of the vector meson 
to the tensor current, /y : 



/^(Ai)=zf (At) {\+bjam. 



q.sea / 



(2.3) 



where e{X)imy 



-(0) 



T^v|V,A) ande(A),-mya/y^'* = {0\a{8^pdpVy - 8vpdpV^)\V ,X). For 



the coefficients bj and cx only perturbative results are available [g]. The QCDSF collaboration has 
calculated the renormalisation constant (which depends on both, a scheme y and a scale pL) non- 
perturbatively in the RI'-MOM scheme [^. We use the 4-loop and 3-loop results for the jS -function 
and anomalous dimension to translate our results into the MS scheme at jJ. =2 GeV [18]]. 
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3. Results 

The results which are obtained for different values of the gauge coupling j3 
and 5.40 can be well parametrised by the following ansatz: 



: 5.20, 5.25, 5.29 



rafx{a, romps,sea, ?"o'"PS,vai) 



roFv(<3) 1 + ao (''0 'MpS.sea)^ + «! (''0 'MpS,val)^ + 



j8o (ronips,sea)'* + i3l "Jps sea"JpS,val + ('"O "IpS,val)'' 



• (3.1) 



Here amps,sea denotes the mass of the pseudo-scalar meson for JCsea = KVai, while for awips.vai the 
valence and sea quark mass might be different. Given the small number of sea quark masses, we 
fixed the value of j3o = 0. In case of /y we find the valence quark mass dependency to be linear in 
(fl;mps,vai)^ and therefore we also fix ^2 = 0. The results of our fits to this ansatz are shown in Fig. ^, 
|and|. As can been seen from these figures, discretisation effects seem to be small, in particular 
compared to the statistical en^ors. For /ps and /y we therefore fit the results in the chiral limit to a 
constant. Only in case of l//v we find much better agreement when fitting to a polynomial linear 
in {a/rof-. 

For j8 = 5.29 and K"sea = 0. 13550, 0. 13590 we determine the decay constants for three different 
volumes x T = 12 x 32, 16 x 32 and 24 x 48, which corresponds to L 1.0, 1.3,2.0 fm. We find 
the pseudo-scalar decay constant /ps to be shifted by 2 — 6% when increasing the spatial extent from 
1.3 to 2.0 fm. This seems to be larger than one would expect from the results of chiral perturbation 
theory [^. No significant finite size effects have been observed for the vector couplings /v and /y . 

In the chiral limit, our result for the pseudo-scalar decay constant is 131(7) MeV, where only 
the statistical error has been taken into account. An extrapolation to the physical pion mass gives 
a slightly larger value of 133(7) MeV, which still agrees very well with the experimental result 
/^"P = 130.7 MeV. Chiral logarithms predicted by chiral perturbation theory, which have not 




0.035 



Figure 2: The upper plot displays the result for ro/ps at j3 = 5.29. The lines show the result of a fit to 
Eq. (3.1). Solid lines are used for the valence quark mass dependency, a dotted for the sea quark mass 
dependency. The lower plot shows ro/ps as a function of {a/ro)- together with a fit to a constant. Both ro/a 
and a/ps have been extrapolated to the chiral limit. 
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been taken into account here, would shift this result down. Extrapolating our results to the kaon 
mass gives us an estimate of /k ~ 155 MeV. This again agrees well with the experimental result 
/k"'' ~ 159.8 MeV. It should be noted however that our estimate for /k is based on simulations 
with degenerate quark masses. For the rho we find l//p = 0.333(12). This is slightly larger 
than the experimental estimate l/Zp"'' ~ 0.27 obtained from the branching ratio BR(t v^). 
Extrapolating our results for the tensor coupling gives /y = 183(4) MeV or, after conversion to 
the MS-scheme, fy^^{2 GeV) = 168(3) MeV. This result is larger compared to previous results 
from QCD sum rules [ p^ ] and from lattice simulations employing the quenched approximation [|l], 
ni], which are in the range of 140 - 150 MeV. 
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Figure 4: Similar to Fig. 2, but for fy^'^^. 
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